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Abstract 

All microscopic correlation functions of the spectrum of the Hermitian Wilson Dirac operator 
with any number of flavors with equal masses are computed. In particular, we give explicit results 
for the spectral density in the physical case with two light quark flavors. The results include the 
leading effect in the discretization error and are given for fixed index of the Wilson Dirac operator. 
They have been obtained starting from chiral Lagrangians for the generating function of the Dirac 
spectrum. Microscopic correlation functions of the real eigenvalues of the Wilson Dirac operator 
are computed following the same approach. 
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I. INTRODUCTION 



The deep chiral limit of QCD with two quark flavors and the intimately related nature 
of spontaneous chiral symmetry breaking is of direct phenomenological interest. Also for 
studies beyond the standard model such as in QCD like theories with many light flavors 
or where the fermions are outside the fundamental representation, the deep chiral limit is 
central. By a remarkable series of numerical and analytic developments it is now possible 
to access the chiral limit by means of lattice QCD. The work presented in this paper is an 
attempt to facilitate the next step to the deep chiral limit by offering an exact analytic un- 
derstanding of the average behavior of the small eigenvalues of the Wilson Dirac operator at 
nonzero lattice spacing, a. The behavior of these eigenvalues is essential for chiral symmetry 
breaking 1-3| as well as for the stability of lattice QCD computations {4]. 

We consider the eigenvalue density of the Wilson Dirac operator in the microscopic scaling 
limit {g-QI where the product of the eigenvalues and the four- volume, V ^ as well as the 
product a?V are kept flxed. This part of spectrum is uniquely determined by js] global 
symmetries, their breaking and the 75-Hermiticity of the Wilson Dirac operator 



(1) 



Because of this Hermiticity relation, the eigenvalues of the Hermitian Wilson Dirac operator. 



-D5 = 75 (-Dm/ + m) 



(2) 



are real. In addition to correlations of these eigenvalues, we will also analyze the real 
eigenvalues of D^r in the microscopic limit. 



In a recent letter [5| and a longer follow-up [10| we have shown how the quenched micro- 
scopic Wilson Dirac spectrum can be obtained from the chiral Lagrangian including order 
a^-effects for Wilson fermions. Although the supersymmetric method used in can be 

applied to any number of flavors, the proliferation of terms makes the method only practical 
for use in the quenched case. Already for one dynamical flavor it becomes rather tedious to 
deal with analytically 11 1. 

In this paper we follow a different path, the graded eigenvalue method, that results in 
simple expressions for any number of flavors with equal quark mass. It is also possible to 
write down compact expressions for all spectral correlation functions. The graded eigenvalue 
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method is based on the observation that the order terms in the (graded) chiral Lagrangian 
can be hnearized at the expense of extra Gaussian integrations. This results in compact 
expressions for microscopic Wilson Dirac spectra for any number of flavors and all correlation 
functions. The method was originally developed to describe transitions between different 



universality classes of Random Matrix Theories [12|-|15|. The result obtained this way is an 
expression in terms of diffusion in superspace where o? plays the role of time. 

All results in this paper are given for fixed index of the Wilson Dirac operator, defined 
for a given gauge field configuration by 



V 

k 



^sign((fc|75|fc)). (3) 



Here, \k) denotes the fc'th eigenstate of the Wilson Dirac operator. The microscopic eigen- 
value density for fixed v gives detailed information on the effect of a nonzero lattice spacing 
on the would be topological zero modes at zero lattice spacing. 

The study of the spectrum also casts new light [s], [lo| on the additional low energy 



constants of the chiral 
theory as developed in 
spacings can be found in 



agrangian which is the backbone of Wilson chiral perturbation 



reviews of effective field theory methods at finite lattice 
24|). By a match of the two- flavor results presented in this 
paper to the microscopic spectrum of the Wilson Dirac operator on the lattice, the value of 
the low energy constants can be measured. The spectrum of the Hermitian Wilson Dirac 



operator in the p-regime of Wilson chiral Perturbation Theory has been discussed in 



25| 



and the results at next to leading order have been fitted to lattice data in j26|. 

The paper is organized as follows. Starting from a chiral Lagrangian for spectra of the 
Hermitian Wilson Dirac operator at nonzero lattice spacing, we derive compact expressions 
for all spectral correlation functions for any number of flavors. In the second part of this 
paper, we obtain expressions for the distribution of the chiralities over the real eigenvalues 
of the Wilson Dirac operator. Some technical details involving Efetov-Wegner terms are 
discussed in Appendix A, and in Appendix B we give explicit expressions for partition 
functions in terms of an integral over a diffusion kernel. 
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II. SPECTRAL PROPERTIES OF THE HERMITIAN WILSON DIRAC OPERA- 
TOR 



The generating function for p-point spectral correlation functions of the eigenvalues of 
the Hermitian Wilson Dirac operator for QCD with Nf dynamical quarks in the sector of 
gauge field configurations with index u is given by (recall that = '-/^{Dw + m)) 



yv 



^ ^ det(D5 + Zk' - ^eklr,) I ^ ' 



The average is over gauge field configurations with index v weighted by the Yang-Mills 
action. For = this is just the Nf fiavor partition function. We will evaluate this 
generating function in the microscopic limit where V oo with 

mV, ZkV, zu% aV (5) 

kept fixed. The axial masses z^ are required when we apply the graded method to obtain p- 
point eigenvalue correlation functions of the Hermitian Wilson Dirac operator. For example, 
from the graded generating function Z'^^j^^-^{m, z, z'\ a) we can obtain the spectral resolvent 

^^,+111(2, m; a) = lim -^Z^ .^^^{m, z, z'; a), (6) 
and the density of eigenvalues, p'^{X^, m; a), of follows from 

P5(A',rn;a) = /$^5(A^, - X')\ = ilm[G5^^,^,|,(-A^)],^o. (7) 

\ k I N^. 

To derive expressions for the correlation functions in the microscopic limit we will rely on 
the graded eigenvalue method. Before deriving the general result, we will first consider the 
case p = 0, which is just the A^/-fiavor partition function. 



The chiral Lagrangian for Wilson chiral perturbation theory to O(a^) was derived in 



16- 



181]. In [5| we obtained the microscopic partition function for fixed index u also to order a^. 



by decomposing the partition function according to 

00 

ZM,{m,e-a)= e"''Z'i,^Xm-a). (8) 

In the microscopic domain, the partition function reduces to a unitary matrix integral 

Z'^.{m,z;a) = [ dU defU e^^^\ (9) 

■' Ju{Nf) 
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where the action S[U] for degenerate quark masses is given by 

777 Z 

S = —J:VTiiU + U^) + -SVTr(f/ - U^) (10) 



-a' 



In addition to the chiral condensate, E, the action also contains the low energy constants 
Wq, W-j and as parameters. Since the and Wj terms can be eliminated at the expense 



of an extra integration [10|, we will only consider the term in the remainder of this paper. 
For the reasons discussed in section VII of [l^ we consider only > 0. To simplify our 
notation below, we will absorb the factor VW^ into and the factor VL into m, Zk and 

aVPF8^a^ mVT.^m, ZkVE ^ Zk and AVS ^ (11) 



Up to a normalization factor, the term proportional to Ws in the action ( ITOl) can be 
rewritten as 



-a2Tr(C/2+C/t2) _ 2Nfa^-a'^Tr{U-U'')^ 



where a is anti-hermitian, and c is a normalization constant. In a diagonal representation 
of a denoted by S* = diag(isi, ■ ■ ■ ,isN^) the partition function with index u is thus given by 
(up to a normalization constant) 



ZL(m;a) = e-^^/-^' / dskA\{sk}) f rft/ det^t/e" ^'^/^^'^'e^T'-^'^^+^^+^T'-^^"^-^) 

J Ju(Nf) 

= j rfsfcA2({sfc})e-^^^'^/^'^'^'det^(m-^)Z]^/{(m2 + 4)V2};a = 0). (13) 
The Vandermonde determinant is defined by 

p 

A(xi, ■ ■ ■ ,Xp) = JJ(xfc - xj), (14) 

k>l 

and an explicit expression for the partition function at a = is given in Eq. (129|) . The 
expression for the partition function will be discussed in more detail for Nf = 1 and 2 in 
Appendix B 



The fermionic partition function has been written as an integral over a diffusion kernel 
times the partition function at a = 0. Next we will show that exactly the same is true for 
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the graded generating functional obtained in j^. The generating function of the microscopic 
p-point spectral correlation functions of the Hermitian Wilson Dirac operator with index v 
is given by 

^]^,+pb(-^,^;a) = yrft/Sdet(z[/)'^eiT^s(-M[c/-c/-])+i^g(2[c/+t/-])+a2Trg(c/2+c/-)^ (15) 
where M. = diag(mi . . . mjV/+2p) and Z = diag(2;i . . . ZNf+2p), and the integration is over 



Gl{Nf +p\p)/U{p), see [271]. We use the convention that TrgA = Tr[A/] — Tr[Ab], with Af 
the fermion-fermion block of A, and Af, its boson-boson block. The definition of Sdet follows 
form the relation Sdet A = exp[Trg logA]. Notice that in comparison to (Q the integration 
over U has been rotated by i so that the convergence of the bosonic integrals is assured 
[5], [l^. A similar rotation is necessary when computing the spectrum of the Dirac operator 



at nonzero chemical potential [28|]. The common origin is the non-Hermiticity of the Dirac 
operator. The manipulations from ([9]) to (fT3l) can be repeated for the graded partition 
function. We start from the identity 



where a is an (A'^^ + p\p) graded "Hermitian" matrix (see fl2T]) ) and c is an integration 
constant. After shifting integration variables a ^ a — Z we obtain (the normalization 
constants will be fixed at the end of the calculation) 

^]^,+p|p(-M,^;a) = e-^^/'^Y da y^t/ Sdet^(zf/)e^^*^('^-^)'/i6'^'+^^^s('^+^)^+^^'-s('^-^)^"\ 

(17) 

We will evaluate this partition function for the {p + Nf \p) graded diagonal matrix 

Z = diag(e, ■ ■ ■ ,e,zi,- ■ ■ ,Zp,z[,- ■ ■ ,Zp), (18) 

and it is understood that the limit e — t- is taken at the end of the calculation. 
The partition function (I9l) and the generating function (|T5l) satisfy the relation 



Z^J(m, 2;;a) = Z'^^{m, -z; a), (19) 



^^;+,|p(-M, 2; a) = Z'-^^XM, -Z- a). (20) 
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For this reason we only consider the case u > below. For z = {Z = 0) the generating 
function does not depend on the sign of u. 
The graded matrix a has the structure 



laj a 
(3 0-5 



(21) 



with Hermitian {p + Nf) x {p + Nf) matrix aj- = aj, Hermitian p x p matrix al = (Tf,, and 
a and /3 are matrices in the Grassmann algebra. The matrix a can be diagonalized by a 



super-unitary transformation 29 1 



, is , 

a = uSu-\ with S = \ I , (22) 

where s = diag(si, . . . , Sp^^f) and t = diag(ti, . . . ,tp) contain the real eigenvalues Sk and t^. 

Next we transform to the eigenvalues of a and the super-unitary matrix U as integration 
variables. The measure is given by 

Nf+p p 

da = 5^^.+p|p(5') Yl dsk Yl dtk du. (23) 

fc=l k=l 

with Berezinian 

rNf+p 



nfc> ; {'^Sk-isi)Yll>i{tk - ti) 

The measure du is the superinvariant Haar measure. 

For degenerate quarks masses, Aik = fn, the integral over u can be performed by a graded 



generalization 
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30( 1 of the Itzykson-Zuber formula, 



° (27rr)<"/+2ri/2 B„,+,|,(X)B„,+H.(-R) ' ' 

Here, X and R are the diagonal representation of ^ and p, in this order. Bq contains the 
contributions due to the boundary terms. They result from the product of the infinity due 
the singularities in the measure and the vanishing result due to the Grassmann integration 
after changing to eigenvalues as integration variables. These contributions can be worked out 



by expanding the eigenvalues of a in powers of the nilpotent terms (see Appendix A). They 
do not contribute to the spectral correlators discussed below and will be further analyzed in 
a future publication. 



Ignoring the contribution from Bq we find the generating function 



Nf+p\p{ 



The integral over U can be expressed in terms of the generating function for a = 0. A 
convenient expression is obtained by observing that the masses Sk + rn and Sk — m can be 
replaced by the same mass i^rn? — S\ for both chiralities. A further simplification results 
from the symmetry of the integrand in the and the ifc variables so that all terms in 
the Laplace expansion of the determinants from the Itzykson-Zuber integral give the same 
contribution. In terms of the a = partition function we thus obtain 



nL=i(4-^ore 



P ^1 Nf „-2Nfa^ 



^^Tr,iis-zf,nea^ f TTT^^l ^ ^n.^pIp ({(4 + m^Y^ {{m' " t^^; a = O) , 



(27) 



where Vandermonde determinant, A{S/8a^) = A(?si/8a^, . . . , isTv^/Sa^), and the prefactors 
result from the limit e — )■ 0. 

We have succeeded in rewriting the a dependence of the generating function as an integral 



over the product of a diffusion kernel and the a = generating function given by 31 



32| 



^Nf+p\pi^i^ • • • , XNf+2p', a = 0) (28) 
llfc=Ar^+p+i^fc \ det[(a;fc)'"%+i_i(xfc)] 



nfii'^ J ^(^1' ■ ■ ■ ' ^Nf+p)^i^'Nf+p+i^ ■ ■ ■ 



Nf+2pJ 



with 



Xg(xfc) = Ig{xk), k = !,■■■ ,Nf+p, 

Igixk) = i-iyK.ixk), k = Nj+p + l,--- ,Nj + 2p. (29) 

We added a tilde to Z'^^_^_p^p because the mass-factors due to the zero modes have been 
amputated which is not the case for the partition function at a 7^ 0. 

The p point correlation function is obtained by differentiating with respect to the Zk and 
putting = Zk afterwards. Only if all factors in the product Yli^'k ~ ^k) are differentiated 



8 



do we get a nonzero result. The remaining factors from the Berezinian B]y^^p\p{Z) canceL 
We thus find the correlator of p resolvents 

g~2Nfa^ 1 f 

Gn,+p\p{zu--- ,zp;a) = (^-a) (i6^^2)(jv,+2p)/2 J dsdtBN^^pip{S)A{S/8a^) (30) 



^2 



This result is universal in the sense that it is completely determined by the symmetries of 
the QCD partition function. Below we discuss explicit results for the microscopic spectral 
density for Nf = 0, 1 and 2. We also check that the partition function for Nf flavors reduces 
to previously derived results for Nf = and Nf = 1. 

A. Explicit Results for the Spectrum of 

In this section we discuss explicit results for the microscopic spectral density of for 
the quenched case and one and two dynamical flavors. 

1. The quenched case 

In this case the Berezinian is given by 

B,,(S) = ^ (31) 
and the generating function reduces to 

z' - Z f dsdt _r/„_L,•,^2_L^+_,'^2l/1fl„2 Cm - is'" 



' IQa'^n J t — IS [m — ty 



Here, 



xZ[|i(Vm2 + s2, -t2; a = 0). 



Zi\i{x,y; a = 0) = ^[yK^+i{y)h{x) + xK^{y)h+i{x)]. (33) 



The resolvent is given by 

Gi\i{z) = Yz Zl\i{^,z,z'-a). (34) 
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Only the term where the pref actor z' — z is differentiated contributes to the resolvent. In 
the microscopic limit, this results in 



Grii(^) 



-e 



16a^7r J t — is [m — t — z) 

xZl\^{^m? - {is + zf, ^Jm? - {t + zf] a = 0). 



(35) 



Note that we also shifted the integration variables s and t by —iz and z', respectively. The 



effect of this shift is discussed in Appendix A 



The quenched microscopic eigenvalue density of follows from the imaginary part of 
the resolvent, cf. Eq. ([7j). We have checked numerically that this result coincides with the 
result obtained from a standard supersymmetric computation in fl. See flliol for plots of 
the quenched density. 

In the a — )■ limit at fixed m and z, the Gaussian integrals in Eq. (133|) become (5-functions 
which can be integrated resulting in 

[z — m) 



lim Z'(,,(m, z, z'\ a) 



Z'(JVm^ - z\ ^/w? - z'^- a = 0). 



[z' — m) 

The prefactor gives the contribution from the zero modes to the resolvent 



z — m 

whereas the second factor gives the contribution of the nonzero modes for a = 0. 



(36) 



(37) 



B. One Flavor 



For Nf = 1 the generating function is given by 



^2|1 



,Z — Z 



64a37r3/2 



^ f dtdsids2{iS2 - ^■si) ^_[^2^(s^+j^)2^(t_^/)2|/^g^2 

z J {t - isi){t - is2) 



X- 



Im — isi) [m — IS2) 



(38) 



{m-ty 



-~ZU{s\ + m'fl\ (si + mY', im' - ^f'^- a = 0), 



where 



Z^m(xi,X2,X3) 



Xn 



X]^d,2\X2 -''1, 



■det 



h{xi) xih+i{xi) xll^^2{xi) 

h{x2) X2ly+l{x2) Xlly+2{X2) 

\ {-IYK,{X^) (-1)-+1X3K.+1(X3) {-irxlK,+2{x^) ) 



(39) 
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To obtain the resolvent we differentiate with respect to z' and put z' = z after differ- 
entiation. An additional minus sign arises because z' is the bosonic source term. We thus 
find 



-2a-' 



dsids2dt{is2 — is 



QAaH^I'^Z'({m- a) J {t - isi){t - iss) 



iZe-[*?+(^2+iz)2+{t-^)2]/16a2 



(40) 



X ill^-^filll^^-^fHl^Z^^ ((,2 + ^2)1/2^ (^2 ^ ^2)1/2^ (^2 _ ^2)1/2. ^ ^ 



(m - ty 

where the S2 and t integration contours are shifted such that S2 + iz and t — z run over 
the real axis. The resolvent is normalized with respect to the one flavor partition function 
given in Eq. ([22]) • As for the quenched case we have checked numerically that the one flavor 
microscopic density, which follows from the imaginary part of Ggi^, is identical to the result 



ll| obtained from chiral perturbation theory using the standard supersymmetric method. 



We refer to 



ll| for plots. 



The small a limit for fixed m and z can be obtained by first shifting S2 — S2 — iz and 
t — > t + z and then expand the nonexpential factors in Eq. (15^ . The integration measure 
can be expanded as 

(is2 - isi + z) 



■ N/ • . =\^ + -{is2-t)] — ^— 
[t + z — lSi)[t — IS2) z t — IS2 



+ 



(41) 



For small a the partition function Z2\i{{s\+m?y/'^ , ((s2 — izY + m^)^/^, (m^ — (t+2;)^)^/^; a 
0) (denoted by Z2\i below) can be expanded to first order in Sk and t 

d ^ 



m "Z'^im; a = 0) + Si ^^211 

dsi 



n=s2=t=0 



S2 -, — ^211 

dS2 



3l=S2=t=0 



■ (42) 



=S2=t=0 



The term linear in Si does not contribute to leading order in a and the remaining terms can 
be written as 

d ^ 



m ''Z'{{m; a = 0) + {t - is2) ^^2|i 



=t=o 



m " Zlim; a = 0) + (t — 182) 



d 

dl' 



Z2\i{m, y/m'^ - z'^, 



.'2^ 



ISi t — IS2 
— V h V 

m m — z 



+ ■ ■ 



The linear term in si in the expansion of the prefactor 

{m — isiY {m — z — is2Y m'^{m — z] 
{m — z — ty {fn — z'Y 

also vanishes after integration. Combining the contributions from Eqs. (Hi 
the small a limit 

Ciz) = + G\m, z; a = 0) + -, 

z — m z 



(43) 

(44) 
we obtain 

(45) 
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where G^{m^ 2;; a = 0) is the resolvent of the nonzero eigenvalues. The additional 1/ z term 
will be canceled by Efetov-Wegner terms which contribute to the real part of the resolvent. 

C. Two Flavors 



In this section we write out the Nf = 2 generating function with dynamical quark masses 
m and index u. With the Berezinian given by 

{is2 - isi){is3 - isi)(is3 - 182) 



we find the generating function 



^311 



z'^z'-z 



{t - isi){t - is2){t - is-i) 

dsids2ds3dt{is2 — isiY^is^ — isi){is3 — 182) 



7r^(16a' 



2^3 



z^ 



(t - isi){t - is2){t - iss) 



y^^-[sl+sl + (s3+izf + (t-z'f]/16a 



2 [m — tsi^ 



[m 



is2Y{m - issY 



2N1/2 



m 



2N1/2 



The partition function for a = is given by 



(m - t)" 

2\l/2 I 2 

m ] 'Am 



1 



t 



2^1/2., 



0). 



X det 



'Zl\\{X\,X2,X-i,X4\a - 2- , 2 _ 2W™2 _ „2W™2 _ ™2N 

( /^.(Xi) x\ly^2{X\) xf/^+3(Xi) 

ly{x2) X2lu+l{x2) Xpu+2{X2) xllu+3{x2) 

Iy{x-i) X3ly+i{x3) xlly+2{xz) xlly+-i{x:i) 

\ {-iyK,{x4) x4(-i)^+i/c+i(x4) xi{-iy^^K,+2{x4) xi{-iy^^K,+;{x4) ) 



(46) 



(47) 



(48) 



The resolvent is obtained by differentiating the factor z' — z with respect to z dX z' = z. 
This results in 



Gzu{z]a) = / dsids2dssdt- 



11^'^ 7r2(16a2)3Z];^.^2("^;«)i {t-ts,)it-iS2){t-iss) 

(m - ty 

XZ3^|,((S? + m^)V2, (,2 + ^2)1/2^ (^2 ^ ^2)1/2^ (^2 _ ^2)1/2. ^ ^ g)^ (49) 

where the S3 and t integration contours are shifted such that S3 + iz and t — z run over 
the real axis. The resolvent has been normalized with respect to the two-flavor partition 
function defined in Eq. (11021) . 
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T ^ 1 I , ^ , I , ^ , I r 




FIG. 1: Spectral density of D5 for Nj = 2, v = 0, m = 3 for a = 0.25 (black curve) and a = (red 
curve) . 




X 



FIG. 2: The spectral density of the eigenvalues of -D5 in the sector with index z/ = plotted for 
Nf = 0, 1 and 2. The value of a = 0.25 and m = 3. The increasing repulsion from the origin for 
larger Nj is clearly visible. 

The microscopic eigenvalue density of with two light flavors of mass m is then given 

by 

p^•^^='(A^m;a) = -lm[G^,^,i-\%^,. (50) 

In Fig. [T]we show the two-flavor microscopic spectral density of as a function of for 
u = and m = 3 and compare the result for a = 0.25 and a = 0. The area below the 
two curves is the same within our numerical accuracy. In Fig. |2]we compare the two-flavor 
result to the one-flavor result and the quenched result with the same parameters. 
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It 
J! 




FIG. 3: The effect of the index on the spectral density of the Hermitian Wilson Dirac operator for 
two flavors. Results for index u = 1 (left) and v = 2 (right) are shown for m = 3 and a = 0.25 (red 
curve) and a = (black curve). The vertical black line marks the position of the v fold (5- function 
due to the exact topological modes at a = 0. Note that the primary effect of a for small a is to 
smear out the (5-function. 

The effect of non-zero index v for A'^^ = 2 is displayed in Fig. [31 Note that the two 
flavor eigenvalue density is positive definite (the square of the Wilson-fermion determinant 
is real and positive) and that the spectral gap cannot close completely on the microscopic 
scale due to the repulsion from the origin (the square of the Wilson-fermion determinant 
vanishes quadratically as an eigenvalue of approaches zero). It would be most interesting 
to compare these analytical predictions to dynamical lattice data, such as those presented 



m 



For small a at fixed m and z we can write 
(zss - isi + z){is3 - is2 + z) 



(51) 



(t - isi + z){t - is2 + z){t - iss) 

The constant term contributes to the real part of the resolvent, and as in the one-flavor 
case, we expect that it will be canceled by contributions from the Efetov-Wegner terms. 
The a — )■ limit of the Nf = 2 resolvent is obtained by expanding the pre-exponential 
factors in Eq. (149|) as in the one-flavor case. 
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III. THE DISTRIBUTION OF REAL MODES 



We now consider the eigenvalues of the usual Wilson Dirac operator Dw- For small 
nonzero values of the lattice spacing a, the eigenvalues, A^, of spread into a narrow 
band around the imaginary axis of the complex eigenvalue plane. The eigenvalues in the 



complex plane make up complex conjugate pairs A'^, (A^)* or are exactly real. See [33| for 
a derivation of these properties of the Wilson Dirac operator. In this section we analyze the 
microscopic spectral correlation functions for the real eigenvalues of Dw- 

The generating function for the p-point correlation function with Nf dynamical quarks 
in the sector of gauge field configurations with index u takes the form 



A+N,fn , ^rr det{Dw + mk) \ , . 

det f [Dw + mf)\ \ — — ■ ) . (52) 



The spectral resolvent for the one point function is 

+i|i(m,m/;a) = lim -^Z](r,+i|i(m/, m, m'; a). (53) 

' m'—>-m am ' 

To be precise, the one point function that corresponds to this resolvent is the distribution 
of the chiralities, sign((A;|75|A;)), over the real modes, A^ G M, 

p^(A^m;;a) = ( 5(Ar + A^) sign((fc|75|A:))\ . (54) 

It can be obtained from the discontinuity of the spectral resolvent across the real axis (see 
section II of [lol ]) 

p^(A'^,m/;a) = ^lm^''j^^^^^{mf,m = \^-ie]a)]^^Q. (55) 

The p-point spectral resolvent is given by 

^jVf+p|p("^i> • • • , mp, m/; a) = lim . . . lim (56) 

xZ](r^,+p|p(m/, mi, ... , mp, m[, . . . , m^; a). 

As in the case of the one-point function, the discontinuities across the real axis give the 
p-point density correlation functions. 

As was discussed in j^l , the generating function for the correlation functions ( 156|) is given 

by 

Z]^^.+p|p(>l;a) = ydf/Sdet(z[/)'^e^^^^(^[^-^"^l)+'^'^^s(c/^+f/-)^ (57) 
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which is just the generating function (1151) for Z = 0. The mass matrix corresponding to 
fl52|) is given by the {p + Nf \p) graded diagonal matrix 

= diag(mj, ■ • ■ ,mf, mi, • • ■ , rUp, m'^, ■ ■ ■ , m'p). (58) 

The first Nf entries are the physical masses and need not be identical. 

In order to derive the p-point function the we start with the identity (instead of the 
identity (|T6|)) 



= ce^^f^" J rfae^'-s'^'/i^'^'+tTrg-lf^-f^-^), (59) 

where a is an {Nf + p\p) graded Hermitian matrix (see Eq. ( 12T|) ) and c a normalization 
constant.. After a shift of a by we obtain 

ZNf+p\piM;a) = e^^^'^Y da y"(it/ Sdet'^(zf/)e^s('^-'^)'/i6a^+tT>s-(^-t^-^). (60) 

The next step is to decompose cr = uSu~^ with S a diagonal graded matrix (see Eq. fl22l) ) 
and perform the integration over m by a supersymmetric generalization of the Itzykson-Zuber 
integral. We find 



Z'(, , (M-a) = / dsdt -MjSf/Sa^ ^ ^ -Mjsf/Sa^ 

,,Tr,[.^H-^^]Aa.^IL|^ ({v^}, {vW};- = O) . (61) 

For degenerate dynamical quarks the above expression can be further simplified. 

d^^e-Mlsl/8a^ ^ A(mi,-- - ,m^^)A(^//8a2,--- ,4^./8a2)e~"^^^'+-^^-/^/"^' 

xdet[e-'"''^'^/S'^'];^ ;^^^,_^i ... _^^_^p+ permutations of (62) 

All permutations of the Sk give the same contribution. For mj, mj^^j^k we obtain 



A(mi,--- ,miv^) -p-p , 



(63) 

The final expression for the generating function with degenerate quark masses is given by 



xe(^/-^)-^[(-^«n|^Z]^,,,, ({ v^}, { vW}; a = O) • (64) 
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In order to obtain nonzero contributions to the spectral resolvent (156|) all rrik in the pre- 
factor have to be differentiated. Below we give the explicit expressions in a couple of cases 
relevant for current lattice simulations. 



A. The quenched case 



The quenched one-point function p^(A^; a) follows from 



a) 



1 



16a^7r / t 



dsdt 



-e 16a^ 



[is + mY 



IS 



{t + my 

xZ'(.^{^{is + m)\ V(t + m)2; a = 0), 



(65) 



after using (155|) . The explicit form of Z^^ at a = is given in Eq. (!33l) . 

The two-point function in the quenched case follows from the discontinuity of 

1 



E2|2(mi,m2;a) 



dsids2dtidt2 



X 



(167ra2)2 

{is2 + m2- isi - mi)(t2 + m2 - h - mi) 
(ti - isi){t2 + 1112- isi - mi){ti + mi - is2 - m2){t2 - 182) 



— ^ [sl+sUtl+tl] iisi + miy{is2 + m2Y 
{ti + miy{t2 + m2y 



(66) 



X ^2^12 ( V(isi +mi)2, y/{is2 + m2y, 



miY, \/ {t2 + m2)2; a 



0) 



across the real mi and m2 axis. Here the a = partition function takes the form 

x^x'^ 1 
Z2|2(a:i,X2,a;3,a;4;a = 0) = 2—— _ ^2v^2 _ ^2n 

/ Iy{xi) Xily+i{xi) xllu+2{xi) xll^+3{xi) 

Iy{x2) X2h+l{x2) Xll^+2ix2) xll^+3{x2) 

{-lyx^ixs) x^i-iy+'K.+iixs) xii-iy^'K,+2{x3) xU-iy^-'K^+six^) 
\ i-iyKix,) x4(-i)^+i/c+i(x4) xi{-iy^^K,+2{x4) xi{-iy^''K,^s{x,) 

The two-point correlation function contains a term due to self-correlations, 



(67) 



X det 



R2ix,y) = (j2 ^\ I, 
This term can be rewritten as 



1 1 \ 1^ 1 1 



k^l 



X + Xky + Xi 



(68) 



Z_-^ ^ _|_ \ , 7, _l_ \, / _ \ .iL^ -r -I- 



x + XkV + Xk y-x\'^x + Xk y + Xk 



(69) 
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which is singular for ?/ — > x if x and y are on opposite sides of the cut of the resolvent. It 



can be shown in general terms 3J] that such singular terms are due to Efetov-Wegner terms 
and are not included in the expression (IMj) . They will be analyzed in a future publication. 
The two-point spectral correlation function 

p^(Ar,Ar;a) = / '^(\^ + ^r)sign((fc|75|fc))5(Ar + Ar)sign((/|75|0)) (70) 

- ( E ^(^r + Af) sign((A;|75|A:))\ / 5{\^ + Af) sign((A;|75|A:)) 

can also be decomposed into sum of self-correlations and genuine two-point correlations 
p;^(Ar, Af; a) = 5{\Y - Af) / 5{\^ + Af)) (71) 



^ 5(Ar + Af) sign((A;|75|A;))5(Ar + Af) sign((/|75|/)) 
J2 5(Ar + Ar)sign((fc|75|A:))\/5^ 5(Ar + Ar)sign((A;|75|A;))y 

k k 

An important observation is that the sign of the chirality drops out in the expression for 
the self-correlations so that the diagonal part of the two-point correlator gives the density 
of real modes 

p;:(Ar, Af = Af; a) = / 5^ 5{X^ + Ar)\ . (72) 



B. One dynamical flavor 

With one dynamical flavor of mass rrif we have 

(rr^ rr^ -n^ = g'"' f dSids2dt{iS2 + m - iSi - TUf) r,2+,2^,2|/ig,2 

z.2|i^m,m^,aj UaH^I^Z-^{mf,a) J {t + m - is^ - mf){t - 182) 

X ^ "^^^ Z^2\^{^l{^s, + mj)\ v/(.i^TH^, « = 0). 

(73) 

Note that the one flavor theory has a sign problem and consequently the one-point function 

p^(A^,m/;a) = ;^Im[S2|i(m - ie, m/; a)],^o (74) 
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changes sign at A'^ = rrif, see also ll| where this function was derived by a direct super- 
symmetry computation. 

The a — )■ hmit atr fixed m andmj is obtained by expanding the pre-exponential factors 
to first order in the and t. This results in 



S211 (m, m/-; a) = — H ^ 

' m m — nif 



S2|i(m,m/; a = 0). + 



(75) 



The l/{m — ruf) term is expected to cancel against the Efetov-Wegner terms. 



C. Two dynamical flavors 

Finally we give the explicit form of the distribution of the chiralities over the real eigen- 
values of Dw in a sector with fixed index u for the physically relevant case of two light 
fiavors with mass rrif. The spectral resolvent can be expressed as 

{is2 — zsi)^(zs3 + m — isi — mf){is^ + m — is2 — rrif) 



X ■ 



[t + m — isi — mf){t + m — is2 — mf){t — is^) 



^^^^.^Is^+sj+sj+t^] jisi + mfy{is2 + mfYjtSs + my 

{t + my ^ ' 

xZ3|i(y^ {isi + m/)2, ^ {is2 + m/)2, ^/ {iss + my, ^/(t + a = 0), 

where the two fiavor partition function in the prefactor is given by Eq. (11021) and the explicit 
form of Zgi^ at a = is given in ( l48l) . 



IV. CONCLUSIONS 



We have obtained analytical expressions for all microscopic spectral correlation functions 
of the Wilson Dirac operator for any number of fiavors with equal quark mass. In particular, 
we have computed the microscopic spectrum of the Hermitian Wilson Dirac operator in the 
physically relevant two fiavor case and the distribution of the chiralities over the real eigen- 
values of the Wilson Dirac operator. The results were obtained from a chiral Lagrangian for 
the generating function of the Wilson-Dirac spectrum using the graded eigenvalue method. 
We have also given expressions for an arbitrary number of fiavors as well as higher order cor- 
relation functions. We have checked that our results for zero and one fiavor are in complete 
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agreement with a previous calculation based on a brute force supersymmetric computation. 
Since these results are based on a chiral Lagrangian that follows from the global symmetries 
of the lattice QCD partition function they can also be derived from a chiral random matrix 
theory for the Wilson Dirac operator with the same global symmetries. This enables us 
to derive additional results using random matrix techniques which we hope to address in a 
future paper. 

The new results give the leading order effect of the lattice discretization on the spectrum 
of the Wilson Dirac operator also in the physically relevant two flavor case. The analytical 
understanding of the smallest eigenvalues of the Wilson Dirac operator can be used to opti- 
mize the choices of parameters in lattice QCD for which the simulation is stable. Our results 
also offer a direct way to measure the low energy constants of Wilson chiral perturbation 
theory. 

Acknowledgments: We would like to thank Thomas Guhr for encouraging us to apply the 
graded eigenvalue method to this problem. We also acknowledge Gemot Akemann, Poul 
Henrik Damgaard and Mario Kieburg for useful discussions and for a critical reading of the 
manuscript. This work was supported by U.S. DOE Grant No. DE-FG-88ER40388 (JV) 
and the Sapere Aude program of The Danish Council for Independent Research (KS). 

Appendix A. EFETOV-WEGNER TERMS 

In this Appendix we illustrate the effect of Efetov-Wegner terms for the Gaussian super- 
integral 



and da is the integral over the matrix elements of a. Clearly, the integral does not depend 
on zi and Z2 so that after a proper normalization of the measure we have 




(77) 



where a and Z are the (1|1) supermatrices 




(78) 



Z{zi,Z2) = 1. 



(79) 
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The supermatrix a can be diagonalized by 



a = u\ ° \u~^ (80) 
t 



with 



1 « 1 /X 

M = exp . (81) 

/3 

We first perform the integral by transforming to the eigenvalues of a as integration variables 
and then perform the integral over m by a supersymmetric generalization of the Itzykson- 
Zuber integral. This results in (Actually this is a special case of Eq. fl26|) where the partition 
function for a = is put equal to unity.) 

levra^ t-is 



Using polar coordinates we obtain 

Z{zi,Z2) = [ rfrci0e*'^+^^e-[^'+^^^'-2^V^-«<^l/i6'^^ (83) 

IQira^ J 

with 

V ^2 - ^1 

The integral over </) is a modified Bessel function 

^(^1, ^2) = ^^^^^ / rfr/,(ry/.i-.f/8a^)e-[^-^+^^^^^l/i6'^l (85) 

Using that 



00 



dxe-""'ji(/3x) = -(e^'/^" - 1) (86) 
we obtain 

Z(zi, Z2) = -e-[^2-^?l/i6»' + 1, (87) 

which disagrees with Eq. fl79p . The missing contributions are the Efetov-Wegner terms 
which arise due to nilpotent terms at the singularity of the measure. Below we will evaluate 
these terms by regularizing the singularity. 
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We regularize the integral ( !77|) over the matrix elements of a by introducing the factor 

^(Va2 + 62_e). (88) 



Writing out the decomposition ( IHOjl we obtain 

a = s — iaf3{is — t), 

b = t + aPiis - t) (89) 

so that 

af3{is — ty 



V^2^=v/i^- -^v- (90) 



and 



^(Va2 + 62 _ e) = e{^/7^Tfl - e) - 5(7^2^ _ ^^^^(^^^_ ^g^^ 

V S2 + ^2 

The measure is given by 



2 



dsdtdadB . . 

da = — ^ (92) 

(t-Zs)2 ^ ^ 



resulting in 



2n J [t — is)^ 

Expanding the nilpotents in the exponent reproduces the result obtained from the Itzykson- 
Zuber integral which does not have to be regularized. We thus find 

IGira^ I t — is 



This is not the end of the story. Because of a discontinuity in the integrand we cannot 
simply shift the integration over s by —izi. There is an additional contribution from the 
discontinuity across the imaginary s-axis. We have that (see Fig. H]) 



dsF{s) = / dsF{s) (95) 

Ci J C2 



with 



/OO PZl PZl 

dsF{s — izi) + i / dyF{—izi + iy — e)~i / dyF{—izi + iy + e). 
00 Jo Jo 

(96) 
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s 






^ 

C2 





FIG. 4: Shifting the s-integration from Ci to C2 gives an additional contribution from the discon- 
tinuity across the imaginary s-axis. 

Applying this to the integral in Eq. we obtain for the contribution of the vertical part 
of the integration contour 

IGira^ J 

The second term cancels against the Efetov-Wegner term. 

The same derivation can be applied to the calculation of the quenched resolvent. The 
conclusion is that if we shift the z and z' dependence from the exponent to the l/(t — is) 
factor, the Efetov-Wegner term is of the form exp(— (2; — 2')^/16a^) which does not contribute 
to the quenched resolvent. 




Appendix B. DIFFUSIVE PARTITION FUNCTION 

The A'^z-flavor fermionic partition function was derived in Section [III from the chiral La- 
grangian. Including the normalization the Nf flavor partition function in the sector with 
index u is given by 

-2Nfa^ r 

Z^^^Xm,Z;a) = ^^^^^.^j,^./, J l[ds,A{{zSk})A{{zs,/8a'}) (98) 

k 

The normalization factor is such that we recover an identity for a — )■ 0. 
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For Nf = 1 we find 
ZNf=i{m, z; a) 
Using the identity 



-2a^ poo 



/ — oo 



(s2 + m2)^/2 ■ 



(99) 



IS — m 



v/2 



^ d9 



g m sin 9+s cos f 



(100) 



IS + m J " ' ' ./-^ 27r 

the integral over s becomes a simple Gaussian integral, and the one-flavor partition function 
can be rewritten as 

de 



Z'^ ^^{m,z;a) = e-'" / ^e""' 6-"''''''-"'°''+'"' '''' ' (101) 

This is indeed the expression for the one flavor partition derived in jsl- 

For Nf = 2 the normalization factor is given by A/" = l/{7r{lQa'^)'^), so that the two flavor 
partition function reduces to 

4a^ roo poo 



^]vf=2("^i:"^2; a) 



dsids2 e 16^^^ ^ ^ (^■Si) (^■52) 



oo ^ — oo 



nil ~ 1^2 
xZ^,{{{zsi)YM{^S2)Y';(^ = 0), 



(102) 



where 



Z2{Xi,X2) 



det 



/^(a;i) xi/,,+i(xi) 

lv{x2) X2lv+\{X2) 



(103) 



It is also instructive to work out the partition function for N ^ = —1. Using the general 
expression (l26l) we obtain 

„2a2 /-oo 



ZNj=^i{m,z; a) 

Using the identity 

t — ie + m 



I — oo 



^^g-(t-.)Vl6a2 (!^!_J!)!^ (_l)-^^((^2 _ ^2)l/2)_ ^-^Q4) 



v/2 



/oo 
-oo 



t — ie — m ^ 

we obtain after performing the Gaussian integration over t and shifting the s-integration by 

■ni 



At 



N^=-i[m, z; a) = e 



2a^ 



^—us^im sinh s+iz cosh s— 4a^ cosh^ s 



^—us^—im sinh s—iz cosh s—2a^ cosh(2s) 



(106) 
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which agrees with the bosonic part of the result obtained in [5| 
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